Theory Of The One Dimensional SU(4) Spin-Orbital Model 
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The one-dimensional spin-orbital model is studied by means of Abelian bosonization. We derive 
the low-energy effective theory which enables us to study small deviations from the SU(4) symmetric 
point. We show that there exists a massless region with algebraically decaying correlation functions, 
~ cos{j^x)x~^^'^ . When entering the massive phase, the system displays an approximate SO(6) 
enlarged symmetry with a dimerization type of ordering consisting in alternating spin and orbital 
singlets. 
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The interest in spin-orbital models stems from the pos- 
sibility of understanding the magnetic structures of tran- 
sition metal compounds 1^. In most of these materials, 

\ in addition to the usual spin degeneracy, the low-lying 
electron states are also characterized by orbital degener- 

' acy. It is thus believed that the unusual magnetic proper- 
ties observed in many of these compounds should be ex- 
plained in terms of simple multi-band Hubbard-like mod- 
els. Very recently, the discovery of new spin-gapped ma- 

■ terials, Na2Ti2Sb20 ||| and Na2V205 g, have attracted 
renewed interest in the spin-orbital models. These mate- 
rials have a quasi- ID structure Q and are modelled by 
a quater-filled two-band Hubbard model which, in the 
limit of strong Coulomb repulsion, is equivalent to two 
interacting Heisenberg models with the Hamiltonian: 

H. — Ji Si ■ Si+i + J2Ti ■ Ti+i 

' i 

; +k(^s,- (f, • f,+i) (1) 

where Si and Ti are spin-1/2 operators representing the 
spin and orbital degrees of freedom at each site i, and 
Ji^2 and K are positive constants. 

The Hamiltonian (|^) is invariant under independent 
SU(2) rotations in the spin (5*) and orbital (T) spaces. 
It can also be recast as a two-leg spin ladder with a four- 
spin interchain coupling. In the limit K <^ J1.2 this 
interaction, which can be generated either by phonons or 
1 (in the doped state) by the Coulomb repulsion between 
the holes, gives rise to a non-Haldane spin-liquid state 
where magnon excitations are incoherent The phys- 
ically relevant question is whether or not this scenario 
can be extended to larger values of K for which (|^) is 
expected to be of experimental relevance. 

As a matter of fact, we already know that this can- 
not be the case. Indeed, the interesting feature of the 
Hamiltonian (|^) is that at Ji = J2 — K/i it is not only 
SU(2)xSU(2) symmetric but actually has an enlarged 



SU(4) symmetry Q. At this special point, the model 
is Bethe-ansatz solvable ||^ and critical with three gap- 
less bosonic modes; in the Conformal Field Theory lan- 
guage, that means that the central charge is c = 3 and, 
as shown by Affleck the critical theory corresponds 
to the SU(4)i Wess-Zumino-Novikov-Witten (WZNW) 
model. Clearly, there should be a qualitative change in 
the physical behavior of (^) when going from small to 
large values of K. ^From the theoretical point of view 
this situation is striking because it implies that one can- 
not go continuously from weak to strong coupling. This 
is a manifestation of Zamolodchikov's c-theorem which 
states that, starting at K — with two decoupled S=l/2 
Heisenberg chains with the total central charge c — 2 
(two gapless bosons), one cannot flow - in the Renor- 
malization Group (RG) sense - towards the SU(4) point 
which has a larger central charge c = 3. Therefore, the 
physics in the neigborhood of the SU(4) point cannot be 
understood in terms of weakly coupled Heisenberg chains, 
and the general strategy employed to tackle spin ladders 
does not apply here: a new effective theory is to be devel- 
oped. It is the purpose of this work to do so. Below we 
present an effective continuum description of the model 
(|l|) at the SU(4) point, based on Abelian bosonization, 
and derive the low-energy expressions for the spin and 
orbital densities. With these results at hand, we then 
investigate the properties of differing phases occuring at 
small deviations from the SU(4) point. 

-Abelian bosonization at the SU(4) point. We start by 
introducing the SU(4) Hubbard model with U > 0: 

nU^Yl ("^ + H.C.) 

iaa 

+ ^ XI "»'^"»6<t' (l-'^afc^^^')- (2) 
iab(T(T 

Here cj^^ creates an electron with the "flavor" (orbital) 
index a = 1, 2 and spin a =t, i, and Uiaa = cl^a^iaa- It 
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will be assumed that the electron band is quarter-filled 
implying that the Fermi momentum kp = 7r/4ao, where 
flo is the lattice spacing. The spin and orbital operators 
are defined as 



Si = 2 X!"^^aa'?a/3C,a/3, 
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Ti — 2 C^iaoJabC-iba (3) 



where a (t) are the Pauli matrices acting in the spin 
(orbital) space. 

The low-energy physics can be described in terms 
of right-moving (Raa) and left-moving (Laa) fermions 
which replace the original lattice fermion Ciaa in the con- 
tinumn limit [x ~ iao): 

~ Raa (x) CXp (ikpx) + Laa {x) CXp {-\kpx) . (4) 



At this point we introduce four chiral bosonic fields 
^aaR,L using Abclian bosonization of Dirac fermions: 

R{L)aa = Kaa{'2^Tao)-^/^exp{±i^/A^<^>aaBiL))■ Thc 

bosonic fields satisfy the commutation 



relation [^aaR, ^ 



ba L- 



Anticommutation be- 



tween the fermions with different spin-channel indices is 
ensured by Klein factors (here Majorana fermions) Kaa- 
It is then suitable to employ a physically transparent ba- 
sis (cf. Ref. H): 

$c = ($1T + *U + *2T + *2i) /2 
$s ($1T - *U + *2T - 3'2i) /2 
= ($1T + *U - *2T - *2i) /2 
= ($11 -$2T + *2i)/2. (5) 

In the new basis, the total charge degree of freedom are 
described by 4>c, while the non-Abelian (spin-orbital) 
degrees of freedom, are faithfully represented by three 
bosonic fields $a (a — s, /, sf). It is now straightforward 
to obtain the continuum limit of the Hubbard Haniilto- 
nian (|^) which exhibits separation between the charge 
and spin-orbital parts of the spectrum. The charge sec- 
tor is described by a Gaussian model for the field $c 
perturbed by an Umklapp term ~ cos ^/WnKc^c gener- 
ated in higher orders of perturbation theory. Though at 
small U the Umklapp term is irrelevant and the charge 
excitations remain gapless, one expects that on increas- 
ing the Coulomb interaction the non-universal parame- 
ter Kc{U) will decrease and eventually reach the critical 
value Kc{Uc) = 1/2 where a Mott transition occurs to an 
insulating phase |^,^ . Though one certainly expects the 
system to be insulating in thc limit U/t ^ oo [|l|, |lo| , |lT| , 
the question whether a commensurability gap, rric, opens 
at a finite value of U is beyond the scope of perturba- 
tion theory. Very recently, Assaraf et al |jl^ using an 
improved Monte Carlo method were able to show that 
there exists a critical value U ^ Uc ^ 2.8t above which 
nic 7^ 0. Assuming U > Uc, in what follows we shall focus 
on the spin-orbital sector described by thc Hamiltonian: 



E 

a=sj,sf 

G3 V 



T 



2„2 



cos V 47r<I>a cos V 47r<I> 



(6) 



a=^b 



where G3 — — ?7ao/2, and 9a — f^aL — '^aR. are the 
fields dual to ^a- The structure of the last term in 
(^ immediately suggests refermionization of thc three 
bosonic fields $a in terms of six real (Majorana) fermions 
r,a=(l,-,6): 

{e 



ie] 



m 



R{L) 



R(L) 



: exp 



(- 



47r$, 



R{L) 



exp ±iV47r$ 



'fR{L) 



R{L) 



: exp 



(- 



47r$ 



sfR{L) 



(7) 



r]i being another Klein factors. In this representation the 
original SU(4) transformations of the complex fermion 
fields appear as SO (6) rotations on the Majorana sextet 
{e}, reflecting the equivalence SU(4) - S0(6). To get a 
better insight in the symmetry properties of our model, 
let us define the spin and orbital triplets: ^ = {^^^1^) 

and £,t — {^1^:^)- Those transform as vectors under 
spin S0(3)s and orbital S0(3)t rotations, respectively. 
In the Majorana representation, the Hamiltonian 
duces to an SO (6) Gross-Neveu (GN) model: 



re- 



a=l 



ildxii) + G3 




(8) 



with Hi — CijCI- Since G3 < 0, we conclude that the 
interaction term in (^ is marginally irrelevant and the 
model flows towards six decoupled massless real fermions. 
Thus, at the fixed point (G3 — 0), the spin-orbital sector 
is described by the S0(6)i (~ SU(4)i) WZNW model 
with the central charge c = 6 • 1/2 = 3. 

To complete our description of the SU(4)-symmetric 
critical point, we present the continuum expressions for 
the effective spin and orbital densities: 



S = JsR. + JsL + exp (i7ra::/2ao) As 
T = JtR + JtL + exp {i'n-x/2aQ)J\ft - 



H.c. + (-!)"■ -n. 



H.c. + (-l)''''"'nt (9) 



x/ao . 
x/ao . 



Here Js,t are the smooth {k ~ 0) parts of these densities, 
while Afs,t and fis^t are the 2kp = 7r/2ao and Akp — tt/oo 
parts. Notice that Eqs. (||) have a more complicated 
structure than that of the spin density in the usual Hub- 
bard model. The emergence of 4ao-oscillations and the 
corresponding complex fields JVs^t is a consequence of the 
band's quarter-filling. The smooth and 2kp contribu- 
tions can be computed directly from Eqs. (||). We find 
that the chiral vector currents 



J, 



sR{L) 



CsRiL) ^tsRiL), JtR(L) — —7^ ^tR{L) ^^tR{L) 
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are in fact SU(2)2 currents, in contrast with a single 
Heisenberg chain where the smooth part of the spin den- 
sity is a sum of SU(2)i vector currents. The fields J\fs,t are 
nonlocal in the Majorana fermions ^s,t- However, as in 
the two-leg ladder problem , they acquire a local form 
when expressed in terms of order and disorder operators 
(Ta and fia of the six critical Ising models associated with 
the six Majorana fermions. The expressions of Afs.t are 
manifestly S0(3)s,t invariant; here we give only their z 
components: A/"/ = A (iAii/i2<73C4cr5(76 -|- cricr2A*3M4M5M6) 
and A/"/ = A (icricr2A*3M40'50-6 + MiM2cr3(J4M5M6), where A 
is a non- universal constant. At the critical point, the or- 
der and disorder operators have scaling dimension 1/8, so 

the 2/cf densities Afs,t have dimension 3/4. Since the S 
and T densities involve fermionic bilinears, it may appear 
surprising to find Akp contributions n^,*. However, noth- 
ing prevents higher harmonics to be generated in inter- 
acting systems. The structure of fls^t can be anticipated 
by symmetry arguments: these fields should be chirally 
invariant and transform as vectors under S0(3)s^t rota- 
tions. These requirements lead to the following simple 
expressions: = iB A S,sL and fit = iB ^tR A S,tL, 
where B is another non-universal constant. The scaling 
dimension of the ns,t fields is 1. 

-Deviations from the SU(4) point. We are now in a 
position to investigate the properties of the model (|l|) at 
small deviations from the SU(4) point. We shall restrict 
consideration to symmetric perturbations, Ji = J2 = 
K/i -\- G, \G\ K, and postpone the study of a more 
general case to a future publication. Using the low-energy 
representation of the spin-orbital densities, one can ex- 
pand around the S0(6)i fixed point to find: 



T-i ^ - iu/2 l^isR ■ d^isR - • dx£,sL 

- iu/2 (^^tR ■ da:S^tR - itL ■ dx^tL 
+ G3 (ki -\- K2 + K3 + K4 + K5 -\- Kof 
-f G (ki + K2 + H&f + (^3 + K4 -t- Ks)^ 



(10) 



The Hamiltonian (|To|) describes two S0(3)-symmetric, 
marginally coupled, spin and orbital GN models. The G- 
term breaks SU(4)~ S0(6) symmetry down to S0(3)s(g) 
SO (3)/. Notice that all interactions are marginal. This 
is the reason why we have also kept the marginally ir- 
relevant (G3) term which is already present at the SU(4) 
point (Eq. (^). The emerging picture is to be opposed to 
the case of two weakly coupled Heisenberg chains where 
the interchain interaction Jj^ gives rise to a strongly rele- 
vant perturbation (of scaling dimension 1) and thus opens 
a spectral gap at arbitrarily small Jj^ ||^ . The RG equa- 
tions for the couplings in ([l^) are easily obtained at the 
one-loop level: 



G = G2_2GG3, G3 =4G3(G + G3). 



(11) 




FIG. 1. Flow diagram for isotropic couplings. 



In the region B, all couplings are irrelevant and a model 
with initial conditions in B flows towards the SO (6)1 fixed 
point. The system is critical and the correlations func- 
tions g{x,T) = {S{x,t) ■ S{0,0)) = (f(x,r) • f(0,0)) 
display a power law behavior. In the long distance 
limit, g{x, t) is dominated by the contibutions at fc = 0, 
k — 2kp and k = Akp: 



— ^-^[{x + iut) ^ + (x — iwr) ^] 

A2cos(— a;)(x2+uV)-3/4 
2ao 

{-ly^^^B^ix^ + u\^)-\ 



(12) 



the leading asymptotics thus being 5^/2. In the regions A 
and C the interaction is relevant and leads to the dyman- 
ical generation of a mass gap. In the far infrared limit, 
all trajectories flow towards the asymptote L: G = — 2G3. 
There the interacting part Hamiltonian ([l0| ) transforms 
to 



Ti-int = G/2 (ki + K2 + Kq - Ks - Ki ~ K^Y 



(13) 



The flow analysis reveals the existence of three different 
regions: A, B and C, shown in Fig.l. 



Upon the transformation ^tR{L) ~^ iCt-R(-L)i the interac- 
tion jl^ ) is easily seen to acquire an SO (6) symmetric 
form. However, the conclusion that the S0(6) symme- 
try is restored in both phases A and G would be incor- 
rect. The scaling portrait in phase C is similar to the 
crossover sector of the Kosterlitz-Thouless phase diagram 
for the U(l)-symmetric Thirring model where an exact 
(Bethe-ansatz) solution |Q confirms restoration of SU(2) 
up to exponentially small corrections. Using arguments 
given recently by Azaria et al. |l^, we therefore expect 
that restoration of SO (6) is a specific feature of phase 
C, while in the massive region A the nature of elemen- 
tary excitations is more complicated reflecting the exis- 
tence of several energy scales. The development of the 
strong-couphng regime in the SO (6) GN model, describ- 
ing phase C, leads to generation of a fermionic mass. As a 
consequence, (ki,2,6) = — ('^3,4,5) 7^ 0, indicating sponta- 
neous breakdown of translational invariance. Indeed, the 
dimerization operators for each chain, As = {—lySi-Si^i 
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and At = (— l)*Ti • Ti+i, express in terms of the energy 
densities of the two S0(3) spin and orbital GN models: 
As ~ Ki + K2 + K6 and Af ^ K3 + K4 + K5. Therefore 
(As) = —(At) — ±Ao, and the system orders in one 
of two, doubly degenerate, ground states with alternat- 
ing spin and orbital singlets, in agreement with the weak 
coupling results ||] . 

Calculating the exact dynamical correlation functions 
in the massive phases is difficult. While hopeless in the 
broken-symmetry phase A, this task could be accom- 
plished in principle in the symmetry-restored phase C 
since the SO (6) GN model is integrable. The full treat- 
ment which takes into account the Z2-degeneracy of the 
ground state and the existence of topological (kink) exci- 
tations in addition to the fundamental fermion will be 
presented elsewhere. However, since the mass of the 
fermion is smaller than twice the kink mass we expect 
that fermions will dominate at sufficiently low energy. 
Their contribution to the correlation functions can be 
estimated by a mean field approach: 

{S{x, t) ■ S{y, 0)) ^ cos(-^a:) cos(-^y) KoiMR) 

„ (_l)^/'^o^2 Kl{MR) 
{f{x, t) ■ f{y, 0)) ~ sin(-^x) sin(-^2;) K^{MR) 



'2ao 2ao 
- [-ly/'^^B'^ Kl{MR) 



(14) 



where R = \/ {x — yY -\- u^t'^ and Kq{MR) is the real 
space propagator of a free massive fermion. We ob- 
serve that on top of an incoherent background at '--^ tt 
(with weight ^ B^), there is a coherent magnon peak at 
k = 11/2 (with weight ^ A^). This is to be contrasted 
with the situation at weak coupling (K <ti J) where only 
incoherent magnons at fc ~ tt exists M . At this point it is 
worth commenting on the status of the non-universal pa- 
rameters A and B that enter in the expressions of the spin 
densities. The numerical results at the SU(4) sym- 

metric point are in good agreement with the expressions 
(p^). In particular, these results have revealed that the 
peak in the static susceptibility at 2kp is much greater 
that the one at Akp, thus suggesting that A ^ B a.t 
the SU(4) point. This has not to be the case when one 
deviates from the SU(4) symmetric point, and the ques- 
tion that naturaly arises is how, as K decreases, will 
one moves from strong to weak coupling regimes. Since 
our solution for large K captures the properties of both 
regimes, it is natural to make the hypothesis that the 
crossover is encoded in the K dependence of the nonuni- 
versal constants A and B. In the simplest scenario, one 
may conjecture that B{K) will increase as K decreases 
while A{K) should decrease. Since A is found to be zero 
at weak coupling, one may further suspect that it will 
vanish for K smaller than a critical value K]j. Such a 
special point where some oscillating component of the 
correlation function vanishes is called a disorder point 

Let us conclude comparing our results with the recent 



DMRG calculations by Pati et al Our result for the 
phase B is in agreement with the numerical data. In the 
phase C these authors find a doubly degenerate ground 
state which form alternating spin and orbital singlets, 
in agreement with our results. However, they conclude 
that the mass gap opens with an exponent 1.5 ± 0.25 
whereas the bosonization approach predicts that the gap 
is exponentially small with the deviation from the SU(4) 
symmetric point. Moreover, they interpret their datas in 
favor of incommensurate correlations in contrast with 
(p^). In the continuum approach, we found no room 
for incommensuration since the parity breaking ( "twist" ) 
term iSfa {x + oq) • Af^ {x) + H.c, which appears upon 
deviating from the critical SU(4) point and which might 
be a potential source of incommensutations | |T^ , turns 
out to be strongly irrelevant (with dimension 3). Our 
result supports another scenario in which the correlation 
functions contain components at A: '--^ 0, tt and 7r/2, with 
amplitudes depending on K. 

We are grateful to Roland Assaraf, Michel Caffarel, 
and Michele Fabrizio for illuminating discussions. 
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